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Abstract 

The explicit analytical expression for the distribution function of parametric 

derivatives of energy levels (" level velocities" ) with respect to a random change 

of scattering potential is derived for the chaotic quantum systems belonging 

to the quasi ID universality class (quantum kicked rotator, "domino" billiard, 

disordered wire, etc.). 
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It is generally accepted nowadays that the problem of a quantum particle moving in 
a random potential first addressed in the context of Anderson localization has much in 
common with such problems in the domain of Quantum Chaos as " quantum kicked rotator" 
|jl|], the ionization of Rydberg atoms by a microwave radiation 0| and chains of quantum 
billiards (the "domino billiard" ||^). This analogy first suggested in proved to be very 
fruitful for understanding of the phenomenon of the so-called "dynamical localization" 

More recently it was found that there exists a convenient mathematical framework - the 
ensemble of Random Banded Matrices (RBM), see and references therein - generalizing 
the classical Gaussian Ensembles of Random Matrices and allowing for a uniform description 
of the typical features common to the abovementioned systems. Investigation of the RBM 
ensemble helped to reveal a number of universal scaling relations characterizing statistics of 
energy levels and eigenfunctions of all these systems m^fl. Another important feature is 
that the stochastic RBM model can be mapped onto a regular field-theoretical model - a so- 
called nonlinear graded a— model - allowing in some cases for an exact analytical treatment 
and so providing one with a powerful tool of research This nonlinear cr— model turns 

out to be identical to that derived earlier by Efetov and Larkin in the course of study of the 
Anderson localization in disordered wires 

All these facts suggest to introduce a notion of a "quasi ID universality class" of disor- 
dered and chaotic systems. All statistical properties of systems belonging to this class are 
dependent on the only scaling parameter: the ratio x = L/C, between the sample length L 
and the localization length ^. The explicit form of the scaling function was derived analyti- 
cally for the so-called inverse participation ratio measuring the extent of eigenfunctions. For 
other quantities characterizing eigenfunction statistics analytical results are available in two 
limiting cases x ^ 1 (x ^ 1) corresponding to the complete localization (delocalization) of 
eigenfunctions [§,0. As to the statistics of energy levels, only heuristic expressions deduced 
from the numerical data were available so far [Q. 

Quite recently an interesting new developement in study of weakly disordered metallic 
systems and their chaotic counterparts has been made in a set of works by the MIT group 
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10| . Developing earlier ideas from the papers the authors of |T0[ studied the energy 



level motion as a function of some external tunable parameter a. Physically the role of such 
a parameter can be played by e.g an external magnetic field, the strength of a scattering 
potential for disordered metal, a form of confining potential for quantum billiards, or any 
other appropriate parameter which the system Hamiltonian is dependent on. A high degree 
of universality in a " level response" of a generic chaotic system to an external perturbation 
has been revealed. It was found that a set of "level velocities" (LV) Vnipi) = dEn/da , with 
index n labeling different energy levels, can be characterized after a proper normalization by 
universal correlation functions whose form is dependent only on symmetries of 

the unperturbed Hamiltonian and those of perturbations. Related quantities characterizing 
energy level response, such as the distribution of "level curvatures" Kn = d'^Vnjdo? and 
that of "avoided crossings" (local minima of adjacent level spacings) were studied in the 



papers [|ll] , p!2[ . Let us also mention an intimate connection between the "level response" 
characteristics and the system conductance if the role of the perturbation is played by the 
Aharonov-Bohm (AB) magnetic flux , see the detailed discussion of the issue in [?]. 

The most of analytic work on the "level response" characteristics done so far have made 
use of the analogy " Quantum Chaotic Systems - Random Matrix Ensembles" which is now a 



commonly accepted principle in the domain of Quantum Chaos |13| . However, a simulation 
of chaotic (disordered) systems by the classical Gaussian Random Matrix Ensembles (or 
equivalent models) precludes effects of Anderson localization from being taken into account. 
On the other hand, these effects should considerably modify all the results when one deals 
with systems belonging to the "quasi ID universality class" introduced above. This indeed 
was found to be the case in the numerical study of the curvature distribution for the 'periodic 
RBM simulating a disordered ring threaded by the AB flux W\ . 



In this letter I derive, for the first time, the distribution of the "level velocities" for the 
systems belonging to the quasi ID universality class in the most interesting limit of infinite 
sample length when the role of the localization effects is expected to be maximal. As the 
particular model for unperturbed system I use the ensemble of nonperiodic Hermitean RBM 



simulating a quasi ID system ( an isolated piece of wire, an irregular billiard chain, etc.) 
subject to a magnetic field. The class of perturbations considered corresponds to a slight 
random change of scattering potential within the wire. For this reason obtained results are 
not valid for systems with periodic geometry and AB flux playing a role of perturbation. 
The latter case requires therefore for a separate investigation. 

Let us consider an unperturbed chaotic or disordered isolated quasi ID system of finite 
size L having N oc L energy levels En, n = 1,2,...,N and described by a Hamiltonian 7i 
whose statistical properties are adequately simulated by those of the RBM ensemble. Let us 
study the "level response" to a perturbation STi. = aV, with a being a small parameter. For 
this purpose it is convenient to introduce the resolvent operator R^{E) = [E—T-C—aV±ie]~^, 
with e being a positive infinitesemal. Then one has the following selfevident identity: 

\im eTrRt(E)TrR-(E) = limV — ^ — (1) 

^-0 °^ ^ ^ ^^o^[E-E40) + ie][E-Em{a)-te] 

It is well known that for the chaotic quasi ID system of any finite size the spectrum 
consists of the set of nondegenerate levels, the probability V{6E) for two adjacent levels to 
be separated by a gap 6E tending to zero when 6E (see, e.g. the so-called "Izrailev 
distribution" and the recent analytical results by Kolokolov on the related subject 0). 
Let us now impose the following requirements: (i) the limiting procedure in eq.(j^) to be 
performed prior to the thermodynamic limit L —>■ oo and (ii) the parameter a oc e ^ 
when performing the limiting procedure. Then the only nonvanishing contribution to the 
double sum in eq.([^) is given by pairs of coinciding level indices n = m. Introducing the 
notations /i = 2e/a and Vn = dEn/da \a=o one obtains: 
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K{^) = \imeTrR^{E)TrR,^^^{E) = ^Y. ^ ^'^T s{E - E^) (2) 

Performing now formally the averaging over the ensemble of Hamiltonians Ti. (denoted 
by the angular brakets (...)) and introducing the mean level density p{E) one immediately 
finds the relation between K{fi) and the distribution V{v) of the "level velocities": 

—Re{K{fi)) = TTfi^ ^ \ dv = 7ifi e-^'^dk V{v) cos kvdv (3) 

lyp J —oo V -\- jJ, Jo J —oo 
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The relations presented above are actually valid for an arbitrary chaotic or disordered 
system. Let us now specify the Hamiltonian 7i as being a N x N Hermitean random 
matrix with independent gaussian distributed entries 7i(z, j) with zero mean and the variance 
{H*{i,j)H{i,j)) = a ~^ Here the function a(r) is assumed to vanish at 

r — s> oo at least exponentially fast, the parameter b (assumed to be large: 6^1) defining 
therefore the effective bandwidth of RBM. 

In order to perform the ensemble average I employ the Efetov's supersymmetric approach 
The detailed exposition of the method as applied to the RBM ensemble can be found in 
and is not repeated here. After performing all the necessary steps the problem is mapped 
onto a ID nonlinear a— model with the action S[Q] = Sq[Q] + 5S'[(5], where 

SolQ] = ^ E Str{Q^ - g.+i)' + tTcpe ^ StrQ.A 

1 / '?iTrnf\^ '"^'^ 1 — A 

5S[Q]=1:T=it[ -] E V(zi,Z2)V(z2,Z3)...V(.,,^i)5tr \{ Q.„-— . 

Here Str stands for the supertrace p|, the 4x4 matrices Qi belong to the graded coset 
space f/(l, l/2)/[/(l, 1) x f/(l, 1) and A = diag(\, 1, -1, -1). 

The nonlinear a— model with the action 5*0 [Q] describes the statistical properties of the 
unperturbed quasi \D system and it was intensively studied in [§[0. The main parameter 
is the coupling constant 7 expressed in terms of the RBM parameters as follows: 7 = 
(vrp)^ X^r oc 6^ It defines the only characteristic spatial length scale due to 

disorder: the localization length ,^ oc 7. On the more formal level it plays the role of the 
correlation length of the matrix field Qi. That means that the matrices Qi and Qj can 
be considered as equal to each other as long as | i — j |^ 7. Let us now make a natural 
assumption that the spatial structure of the perturbation V is of the same type as that of the 
unperturbed Hamiltonian 7i, i.e. V(z, j) vanishes sufficiently fast as long as | z — j |^ 6. We 
will call such a perturbation V the generic one. Then in view of the relation 1 ^ 6 ^ 6^ oc 7 
one can put Qi^ = Qi^ = ... = Qi^ in the expression for 55* [Q]- 

From the previous experience one can anticipate that the main contribution to the 
integrals over the matrix field Qi is coming from the asymptotic domain Qi ~ 1/(67) 



long as e — > 0. As it will be clear afterwards, the essential values of the parameter /i are of 
the order of /i ~ max{'y~^/'^ ^N'^^"^}. For a generic random perturbation one notices that 
Vo = (V^)(z,i) = Y.rn I V(i,m) p is a deterministic (selfaveraging) quantity of the order 
of unity, whereas V(i,i) is a random quantity with zero mean and variance of the order of 
(V(i, lY) ~ Vo/& ^ Vq. Combining all these estimates together one arrives at the final form 
for the effective action of the nonlinear a— model: 
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S[Q] = SolQ] - \ V.StvY^Q^il - A)g,(l - A). 



(5) 



In view of the local- in-space structure of the last term in eq.(^ the corresponding integral 
over the matrices Qi can be performed by the same transfer-matrix method used earlier 
with minor modifications. As the result one finds j^{K[jj)) = ^I{x) where 

y 

with both functions Y^^\T^y) and y*^^-*(r, ri, y) satisfying the same differential equation: 



I{x)= frfn r^dr r^Y'^'\x-T-n;y)Y^'\T,n;y) (6) 
Jo Jo Jo y 



dr ' ^ dy"^ \ Ag 



Here the "scaling" notations g = — 

Tcp 
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1/2 N 

and X = — were introduced for the sake of 
27 



2Vo 

convenience. The equation eq.(0) should be supplied with the "initial" conditions: 

r«(^ = 0;y) = l , Y^^\T = 0,n;y)=yY^'\Tv,y), (8) 

the first one corresponding to the elastic reflection of the quantum particle at the sample 
edges and the second one related to the details of the transfer matrix method 0. 

In order to be able to deal with the function I{x) efficiently it is more convenient to 
consider its Laplace transform Il{p) = e~^^I{x)dx. One finds: 



Il{v)= r -Yj^'\p;y)Yj^'\p;y) (9) 

Jo y 

where the functions Yl^'^'^'^\p,y) satisfy the system of two ordinary differential equations: 



<3A"\p:y) = -\ , QvY'i\v:y) = --Yt\p;y) ; = 1^ - f A + - + 4V (lo) 
r y <9r vV y y 



One can write down the explicit solution to these equations noticing that the two Whittaker 
functions W_g^!i{y/g) = (piiy) and M_g^!i{y/g) = (p^iy) are eigenfunctions of the operator 
Qp corresponding to the eigenvalue = —\{Ap+ 1 — k^) and the corresponding Wronskian 
is equal to Wi^ = —r[K + l]/r[g + {K + l)/2], where r[z] stands for the Euler gamma-function. 
Thus, one finds Y^\p\y) = Tl{yYl^^\p; y)} , where 



y!^'\p; y) = ^r[l + g]W.g,,i2{ylg) + 7^{l - r[l + g]W.g,,i2{ylg)} 



(11) 



and the action of the operator TZ on any function f{y) is defined by ( z/ = y/Ap + 1) : 

If ry A -y roo rl 7 \ 



(12) 



Expressions eqs. ( ]T0| - p!2| ) provide a formal possibility to find the function Il{p) for an 
arbitrary value of p. Actually, however, managable expressions could be extracted only 
in the limiting case p physically corresponding to the system length L being much 
larger than the localization length ^, i.e. x (x L/C, ^ 1- The main simplification occurs if 
one notices that one can neglect the second term in the expression for Y^^''{p;y), eq.(^), 
provided p ^ 0. As the result the eqs.@-(|l^) can be presented in the form: 



„ r2[l+g] .oo dz 

-g p2 Jo z 



c = 



(13) 



hip -> 0) 

£Y('\z) = lW.g,^/,{z) 

Differentiating the identity CW_g^i/2{z) = over the parameter g and using the condition 
Y^'^\z ^ 0) = one obtains: 



1 + ^(1 + .) 



W. 



■3,1/2 



[z) where ip{z) 



91nr[^] 
dz 



(14) 



Substituting this expression into eq.([l3|) and remembering the relation between {K{n)) 
and I{x) one finds: 
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dg 2 dg"^ 



2 io 



k/2 



sinh(A;/2) 



(15) 



1/2 



Introducing now the scaled "level velocity" Vg = f (ivj^j ' comparing the eq.(|T3p 
with the eq.(^ one restores the LV distribution function V{vs) from its Fourier transform: 



^, , dk 
V{vs) = / —coskv, 
Jo Itt 



k/2 



sinh {k/2) 



sinh^ (nv 



TV 

{nVsCoth{nVs) — 1} (16) 



This expression gives the exphcit form of the LV distribution for the case of long quasi ID 
disordered or chaotic system and is the main result of the present Letter. 

Let us briefly mention that in the opposite limiting case of short systems whose length 
L ^ ^ the effects of localization play no role and one easily reproduces from eqs.(|^-|^) the 
Gaussian LV distribution typical for the chaotic systems studied in earlier papers |lT0| , p!2| . 
solution to the eq.(7) in the domain 0<r<x^lis given by the expression: 

Y{t; y) = Y{r = 0; y)exp-T [y + (17) 

that immediately produces the required result when substituted to the eqs.(H,^. 

Comparing the two limiting cases one concludes that the "level velocities" fluctuate much 
stronger when eigenf unctions are localized: (i) the probability to find values of LV exceeding 
the typical value (f^)^/^ decays in the case of extended states like e"'^''^, i.e. much faster than 
a simple exponential typical for localized states, see eq.(|TB|) and (ii) the mean square (f^) 
is proportional to the inverse localization length 1/^ oc I/7 when localization takes place, 
i.e. is much larger than the value of the order of inverse system size 1/L oc typical for 
systems with extended states and the same number of levels N . To this end it is interesting 
to note that in the papers [10| the quantity (v^) was called the "generalized conductance" in 
view of its meaning for the AB case [?]. The results obtained in the present paper suggest 
that the level response of a disordered system subject to a random perturbation is rather 
related to the so-called "inverse participation ratio" which is inversely proportional to the 
eigenfunction extent. A more detailed discussion of this issue will be published elsewhere 

[?]• 

It seems to be interesting to check all this predictions by a direct numerical simulations of 
the systems belonging to the quasi ID universality class. In particular, for a chain of chaotic 
quantum billiards ( the "domino billiard" |^) one should be able to observe a substantial 
increase in LV fluctuations when passing from the regular chain to an irregular one. 
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